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Abstract 

Following the Renormalization Group scheme recently developed by 
Pietronero et al, we introduce a simplifying strategy for the renormalization 
of the relaxation dynamics of sandpile models. In our scheme, five sub-cells 
at a generic scale b form the renormalized cell at the next larger scale. Now 
the fixed point has a unique nonzero dynamical component that allows for a 
great simplification in the computation of the critical exponent z. The values 
obtained are in good agreement with both numerical and theoretical results 
previously reported. 
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The concept of Self- Organized Criticality (SOC) introduced by Bak, Tang and Wiesenfeld 
(BTW) has attracted a wide interest to understand a class of dynamically driven systems 
which self-organize into a statistically stationary state characterized by the lack of any 
typical time or length scale. Numerical results for systems displaying SOC behavior have 
been extensively reported but only a few theoretical approaches are known to be in 

agreement with numerical simulations in all dimensions. The major source of difficulties in 
the study of SOC systems lies in their inherent complexity that makes the models analytically 
tractable only in a few cases. The Abelian version of the BTW sandpile model, early 
addressed by Dhar [Q, turned out to be one of these exceptions. 

Recently, Pietronero, Vespignani, and Zapperi || developed a new type of real space 
renormalization group approach for dynamically driven systems, able to describe the self- 
organized critical state of sandpile models by defining a characterization of the phase space 
in which it is possible the renormalization of the dynamics under repeated change of scale. 
In addition, it is also possible to compute the critical exponents analytically || . The method 
also reveals the nature of the SOC problems and provides a picture about the universality 
classes of different sandpile models. This scheme of renormalization has been recently im- 
proved by considering increasingly complex proliferation paths and extended to forest- 
fire models 0-0- 

In this report, we follow the renormalization procedure of references but using a 
Greek cross-shaped cell instead of a square cell in the renormalization of the relaxation 
dynamics. The critical exponents that characterize the stationary state are then computed 
and they are found to be in good agreement with previous theoretical results and large scale 
numerical simulations both for the BTW and two state model of Manna. We will see that 
the use of this particular choice of cells simplifies the renormalization equations for the BTW 
model. 

In what follows, we will focus on the sandpile critical height models in two dimensions. 
Sandpile models are cellular automaton defined on a lattice where to each site one assigns a 
variable (to which we will refer as energy). We let the system evolves by randomly adding 
units of energy on the system. When the energy of a site reaches a critical value, it relaxes 
releasing its entire energy to the neighboring sites. The affected sites may become unstable 
triggering new toppling events and so on until all sites are again stable. Three different 
classes of sites can be distinguished: (i) those sites for which the addition of a unit of energy 
does not induce relaxation (stable sites), (ii) those sites for which the addition of a unit of 
energy provokes they become unstable (critical sites), and (iii) unstable sites that will relax 
at the next time step. Open boundary conditions allow the energy to leave the system. 

In this formalism, we will denote by p the density of critical sites. These definitions can 
be extended to a generic scale b by considering coarse grained variables. Thus, a cell at scale 
b is considered critical if the addition of a unit of energy SE(b) induces a relaxation of the 
size of the cell, that is, the subrelaxation processes span the cell and transfer energy to some 
neighbors. According to || the relaxation process can lead to four different possibilities at 
coarse grained levels: the energy can be distributed to one, two, three, or four neighbors with 
probabilities Pi, P2, P3, and p± respectively. Of course, it is also possible that in certain cases 
the unstable sites at coarse grained scale do not transfer energy to their nearest neighbors 
as well as to consider different proliferation problems. We, as in [||,|| , will not consider these 
cases JT2|] . Then, the probability distribution is defined by the vector 
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P = (Pl,P2,P3,Pi) (1) 

with the normalization condition J2t=iPi — 1- 

So, the properties of the system are fully characterized by the distribution (p, P) at each 
scale. The relation between p and ~P can be derived by noting that in the stationary state 



the inflow of energy equals the flow of energy out of the system ||13|| . This implies 

1 
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P'^^JiS . (2) 



which allow us to evaluate the stationary distribution of critical sites at each scale k of coarse 
graining. 

Now, we define a renormalization transformation for the relaxation dynamics. We will 
use a cell-to-site transformation on a square lattice, in which each cell at scale b^ k ' is formed 
by five sub-cells at the finer scale (see Fig. 1). We have chosen this type of cells for 

two reasons: one, because it implies the use of greater cells formed by five sub-cells at the 
finer scale, that is, when we scale up, five sub-cells form a new one at the larger scale; and 
second, one is intuitively tempted to follow the geometry of the relaxation that takes place 
in numerical simulations of sandpile models with energy transfer to N, E, S, W neighbors 



15 



The length scaling factor is then = y/5 (see Fig.l). Therefore, at a generic scale U k \ 
each cell is characterized by an index a, ranging from one to five, indicating its number of 
critical sub-cells at the smaller scale fr^ -1 ). In order to ensure the connectivity properties of 
the avalanche in the renormalization procedure, only those configurations with three or more 
sub-cells at scale can span the cell, transferring energy to i neighboring cells. Thus, 

the starting relaxation processes Pi at scale &( fc_1 ) are renormalized in the correspondent 
process p^ at scale . Besides, it has been shown |HJ that site correlations are averaged 



out in the stationary state. Therefore, taking into account this fact and the spanning rule 
we can write the weight of each configuration a in the stationary state as: 

W (a=3) =2p 3 (l-p) 2 

W {a=4) = 4p 4 (l - p) (3) 
W {a=5) = p 5 

Eq. (Q) gives the probability that a cell at scale b^ k ' has the corresponding number of critical 
sub-cells at scale b^^. 

As an example of the general procedure, in Fig. 2 we have drawn a series of relaxation 
processes p\ — > p\ — ► P2 at scale that contributes to the renormalization of p * at the 

larger scale b^ k \ starting from a configuration of a = 3 critical sub-cells. The process consists 
of the following relaxation events that span the cell from left to right satisfying the spanning 
condition. First, the unstable sub-cell on the left relaxes toward the other critical sub-cell 
(the center one, Fig. 2b) which occurs with probability (l/4)p^ fc , where the index (k — 1) 
denotes that the relaxation takes place at scale M fc_1 ). Second, we consider the process 
in which the new unstable sub-cell also relaxes toward the sub-cell on the right through 
another pi process (Fig. 2c) which again happens with a probability {1/A)p < f'~ 1 \ Finally, 
the sub-cell on the right has become unstable and transfers with probability (2/3)p^ ^ 
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two units of energy one inside and one outside the original cell of size (Fig. 2d). The 

(k) 

series of processes described contributes to the renormalization of p\ . Nevertheless, it is 
necessary to note that the relaxations displayed in Fig 2a-2d are not all the processes that 
contribute to the renormalization of Pi through a pi — > pi — > P2 series. Fig.2e shows a 
p 2 relaxation event that, although involves two neighboring sites outside the original cell of 
size b( k \ also contributes to the renormalization of Pi with probability (l/6)p2 • This is 
a new characteristic inherent to the cell-to-site transformation chosen. Now, if we take into 
account all the processes that lead to p[ k \ for a = 3, one gets 

p<?> = 5 { Qp?- 1 ' + \pV ] + P t 1} ) (\p?- 1} ) (Ip?- 1} + Ipt 1 ' + ip?" 1 ') } 

+l{(\pt i) + \p { r ] + -/t i, +pt i} ) QpS*- 1 ') (lp ( r } + lpt i} + ip?- 1 ')} <4) 

where in Eq. (|4]) the factors | and | refer to the multiplicities of the configurations (see 
Fig.3). 

In a similar way (though much more complicated), one obtains expressions for p\ k \ 
(i = 2, 3, 4), for a = 3 and imposes the normalization condition Yh=i Pi — 1- The procedure 
is repeated taking into account the configurations with a = 4 and a = 5 critical sites and the 
renormalized probabilities at level k are finally derived by averaging over the configurations of 
different a- values including their statistical weights W a (p^ k ~ 1 '). Therefore, the probabilities 
Pi at length scale will be given by 

pf ] = jlW a {p^)pt l \a) (5) 

with W a {p^ k - l) ) and 

pi k i) given by Eq. (|3|) and Eq. (|2]), respectively. As the computation of 
the probabilities p\ k ^ in Eq. ([5]) is rather lengthy and cumbersome, we have developed a C- 
code to compute all the polynomial term coefficients that contribute to the renormalization 
transformation. 

Now, we proceed to explore the scale-invariant behavior of the model by finding the 
fixed-point solution p\ k ^ = p\ k \ In order to do this, we start from the shortest length 

scale characterized by (p(°\p^) and study how it evolves under repeated iteration of the 

transformation (§). For the two state model of Manna the parameters (p^°\p^) are 
(p(°> , 0, 1, 0, 0) whereas for the BTW sandpile we have (p^°\ 0, 0, 0, 1). Here, the initial value 
of the density of critical sites p^ is irrelevant for the dynamics since the system evolves to 
a fixed point whatever be the distribution of critical sites at the small scale dynamics. 

As in Refs [[|||, both models have the same fixed point, but here there is an important 
difference in relation to the value of the fixed point. We obtain for the fixed point the 
value (p*,l?*) = (^,0,0,0,1), that is, in the BTW model one starts from the fixed point! 
This is not indeed the case for the two-state model of Manna for which we need to iterate 
Eq. (|5|) more than twenty times to reach the same fixed point. We believe that this is a 
consequence of our renormalization strategy for the relaxation dynamics and constitutes a 
great simplification in the calculation of the dynamical exponent z. In fact, we were ex- 
pecting the existence of a critical fixed-point value different from that reported in references 
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|5||| although the critical exponents should be very close since they are determined by the 
properties of the system at large scales. 

The exponent r that characterizes the power-law avalanche size distribution can be ob- 
tained following the procedure of ||. Consider the probability K b ( k -i) 6 o) that the relaxation 
processes that are active at scale fe^ -1 ) do not extend beyond the larger scale b^ k \ This is 
expressed as || 

fe( fe ) 

/ P ( r ) dr ( h (k) \2(l-r) „ n . 

* = =1"(^4) =l-(V5) 2(l - T) . (6) 

Eq. (||) also satisfies 

K = p*(l - p*) + p*(l - p*) 2 + pg(l - p*) 3 + pl(l - p*) 4 (7) 
Then, the exponent r is given by 

1 llog(l-iT) 



2 log(V5) 



1.235 . (8) 



This value of r is in very good agreement with the value obtained in |5],|6| and with large-scale 
numerical simulations which give r = 1.27 for the two-state model of Manna, and r = 1.29 



for the BTW sandpile model [17 



A second independent critical exponent can also be computed. This is the so-called 
dynamical exponent z that relates the spatial scale r to time scale t through the power law 
t ~ r z . As pointed out in |J, the calculation of z could be an enormous and laborious task 
because the knowledge of the fixed point value is not sufficient and we have to know the 
complete form of the renormalized dynamics. Nevertheless, as we said before, the use of 
our larger cells in the renormalization transformation leads to a fixed point with a unique 
nonzero component in the vector . This constitutes a great simplification in the derivation 
of the complete structure of the renormalized dynamics. In what follows, we will derive at a 
glance the dynamical critical exponent for the BTW sandpile model. In order to obtain the 
dynamical exponent we have to calculate the average number < t > of noncontemporary 
processes at scale fo^™ 1 ) needed to have a relaxation process at the larger scale U k \ which is 
related with z through 

In < t > In < t > 

In Fig. 3 we have depicted the possible starting configurations for the different values of a. 
It is also shown the time steps needed to have a relaxation process at the larger scale. Such 
a simplification in the calculus is possible because we have to consider only the relaxations 
that contribute to the renormalization of P4 at larger scale. As can be seen, we need two time 
steps for the symmetric configurations (those in which the initial unstable site is located at 
the center of the cell) and three for the non-symmetric configurations (those in which the 
initial unstable site is located in one of the critical boundary sites of the cell). Therefore, 
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<t>= EWMV (a)WM (10) 

a 

where t'(a) is the weighted average of time steps taking into account the different additional 
statistical weights due to multiplicities u in each configuration a (see Fig. 3). Now, evaluating 



Eq. (10) at the fixed point we obtain, 

2 = 1.236 (11) 



The value ( [TT]) is in remarkably good agreement with the numerical result z = 1.21 and 
with the exact value 2 = 5/4 ||18|| . The other critical exponents can be derived from scaling 



relations ||19|| . Table I summarizes the values of the critical exponents obtained for the 
BTW sandpile model and those reported by previous renormalization scheme and numerical 
simulations. 

In this report, we have introduced an alternative renormalization strategy that simplifies 
the analytical derivation of the critical exponents that characterize the dynamics of sandpile 
models. By using larger cells, formed by five sub-cells of the finer scale, we obtain a fixed 
point with a unique nonzero dynamical component which allow us to derive the whole form 
of the renormalized dynamics in a more direct and simple way. The values of the exponents 
here obtained are in good agreement with those previously reported. Besides, as in similar 
analytical predictions, the two-state model of Manna and the BTW sandpile model belong 



to the same universality class [^0j. The results confirm the robustness of the renormalization 
group approach. 

It is a pleasure to thank A. Vespignani for stimulating discussions and G. Caldarelli for 
useful correspondence. Y.M would like to thank the AECI for financial support. 
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TABLES 



TABLE I. Values of the critical exponents for the BTW sanpile model (d = 2). We have 
included the values obtained from large scale simulations and those derived in [6]. 



Method 


T 


Q 


A 


z 


RG § 
Simulations [17, || 
This paper 


1.253 

1.29 

1.235 


1.432 

1.38 

1.38 


1.506 

1.44 

1.47 


1.168 

1.21 

1.236 



s 



FIGURES 



FIG. 1. Greek cross-shaped cell used in the renormalization procedure. It is displayed the 
central sub-cell (encircled dot) and its four nearest neighbors (black dots). The length scaling 
factor is 

FIG. 2. A series of relaxation processes p\ — ► p\ — ► p2- Open dots represent stable sites, filled 
dots critical sites, and encircled dots unstable sites. Note that the last relaxation affects only one 
neighbor despite of having two outward arrows (see also Fig.l). 

FIG. 3. Full set of possible initial configurations of critical sites and their multiplicities w. We 
have only depicted the configurations that fulfill the spanning rule (a = 3,4,5). The t's refer to 
the noncontemporary time steps needed to have a relaxation that span the whole cell. The indices 
s and a stand for symmetric and non-symmetric configurations. 
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